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SUMMARY 

With the advent of monolithic microwave Integrated circuit (MMIC) tech- 
nology, the phased array has become a key component In the design of advanced 
antenna systems. Array-fed antennas are used extensively In today's multiple 
beam satellite antennas. In this report, a computer program based on a very 
efficient numerical technique for calculating the radiated power (Romberg 
Integration), directivity and radiation pattern of a phased array Is described. 
The formulation developed Is very general, and takes Into account arbitrary 
element polarization, Eand H-plane element pattern, element location, and 
complex element excitation. For comparison purposes sample cases have been 
presented. Excellent agreement has been obtained for all cases. Also Included 
In appendixes A and B are a user guide and a copy of the computer program. 


INTRODUCTION 

One of the most Important radiation characteristics of an antenna system 
Is the directivity. Accurate determination of this parameter Is essential for 
the analysis and design of advanced antenna systems. In general, for most 
commonly used array elements such as open-ended waveguides, horns, or micro- 
strip patch antennas, an analytical expression of the (cos e^-type can be 
used to properly tailor actual element patterns (ref. 1). Experimentally 
measured element patterns could Include mutual coupling, which may be signifi- 
cant In large arrays. Expressions for the directivity and array radiation 
pattern of a single element and a rectangular array using (cos ©JQ-type 
element patterns have been reported (refs. 2 to 11). 

It Is the purpose of this work to generalize these results and obtain an 
efficient numerical technique for computing the directivity and the antenna 
radiation pattern of the generalized array. The generalized array characteri- 
stics used In this report Includes arbitrary element location, element pattern 
((cos e) c l-type, other analytically descrlbable functions or experimentally 
measured), element polarization and element excitation. 


ARRAY RADIATION PATTERN 

The geometry of the generalized array Is shown In figure 1. Given the 
array geometry and element characteristics, the generalized array problem can 
be defined as: (1) to determine the power radiated and directivity at a given 

observation point (this Is usually taken In the far-fleld zone), (2) to deter- 
mine the co-polarization and cross-polarization component of the electric field 
(using Ludwig's criterion (ref. 12)). In solving this problem, two sets of 
coordinate system are used. Figure 2 shows a typical element coordinate 
system and the reference coordinate system, with the z-axls In the same 
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direction. For an array of H elements located arbitrarily In the reference 
coordinate system (fig. 1), the m th element radiated field Is given by 
equation A1 .1 : 
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for 0 < e 1 < ir/2, where 

I m element complex excitation coefficient 

u Em* u Hm element E and H plane pattern 

a m ,bm,% element polarization parameters (see table I) 

K wave number 2ir/\ 

r^,e',<p' spherical coordinates In the element coordinate system 


(Al.l) 


The element pattern U[r m ,Un m can be described with an analytical 
expression ((cos ©Id-type or other functions) or with experimentally meas- 
ured data (discrete). If measured data are used, the pattern may Include 
mutual coupling effects. The polarization parameters In table I are subject 
to the normalization describe by 


The electric field described by equation (Al.l) Is In the element coordinate 
system. The total electric field due to all M elements Is the superposition 
of the electric field of each element of the array. The total electric field 
Is given by 
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where the vector fields E*("r) and E^CT) are defined ‘In the reference coor- 
dinate system. A transformation of t^(’r'm) (eq. (Al.l)) In the element 
coordinate system Into E^("r) In the reference coordinate system Is described 
next. Figure 3 shows a detailed description of these coordinate systems. The 
transformation of coordinates for this problem only Involves a translation. 

The transformation procedure Is outlined as follows. Knowing the observation 
coordinates (r,o,(p) and element location (x m ,y m ,z m ). the observation 
point In the primed coordinate system Is found by using: 


x = r sin e cos (p 
y = r sin e sin <p 
z = r cos © 



(A1 .4a) 


( A1 .4b) 


2 
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0 ' = cos T 


(A1 .4c) 


9 1 = tan 1 ^r 

Wlth equation (A1.4c) computed, all parameters on equation (Al.l) can be 
calculated. The vector transformation Is obtained by using: 
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(A1.5) 


where I_ Is the Identity matrix (3x3). 

Equation (A1.5) transformed E m (r^) Into E^("r) In the reference coor- 
dinate system. This process Is repeated for each array element. Notice that 
no constraints have been put Into equation (A1.3) regarding the observation 
distance. This expression (eq. (A1.3)) Is valid everywhere except at the 
location of the source Itself. This formulation assumes that each pattern Is 
boreslghted In the +z direction. However, the Identity matrix 1 In equation 
(A1.5) could be replaced by a rotation matrix (Euler matrix of transformation) 
to account for arbitrary pointing. 


The array radiation pattern Is usually divided Into two orthogonal polar- 
izations. Equation (A1.3) may be written: 
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which can be expressed as: 
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t(r) = E e + E ® 

U (p 


(A1.5b) 


The orthogonal components described In equation (A1.5b) are the usual 
spherical components. Another way of dividing the electric field Into two 
orthogonal polarization Is by using Ludwig's definition 3 (ref. 12). The 
following polarization vectors are Introduced: 
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R = e(ae^ cos tp + b sin <pj + <p(-ae''^' |, sin tp + b cos tpj (A1.5c) 

C = e^ae - ^ sin tp - b cos tpj + tp^e - ^ cos tp + b sin tp] (A1.5d) 

The reference polarization and cross polarization expressions of E(r) are: 

Reference polarization of f: E„ = T • R* (A1.5e) 

Cross polarization of T: E £ = if • C* (A1.5f) 

With these expressions equation (A1.5a) can be rewritten as: 

t(r) = E r R + E C C (A1.5g) 

The parameters a, b, and ij; can be obtained from table I. 


POWER RADIATEO 

The total power radiated (time-averaged) by the array Is given by: 
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where 

7 x t = -jwyoH Maxwell equation 

ds a r r2 sin e de dtp differential surface area 

S a sphere of radius r 

In the far-fleld of the array (usually taken at 2 a 2 /\, A Is maximum 
array dimension), the power radiated given In equation (B.l) can be simplified 
to: 


rad 
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(Z q Is free space Impedance) 

Substituting equation (A1.3) Into equation (B.2) gives 
2ir ir/2 / M \ / M \* 


rad 


/ / £ 

0 0 \ m=l 
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In general, the above expression does -not have a closed form solution and 
Is evaluated numerically using a Romberg Integration algorithm (ref. 13). 
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Using far-fleld approximations, equation (A 1 . 1 ) can be simplified as follows: 
(This will restrict the observation distance to be only In the far-fleld of 
the array.) 
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E m (r) In equation (B. 4 ) Is given by: 
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With these far-fleld approximations all position vectors "r m are parallel, 
making (§',9') equal (§,9). No further coordinate transformation Is 
required. Substituting equation (B. 4 ) Into equation (B. 2 ) produces: 
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the equation (B.6) can be expressed In the matrix form: 
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The coefficient R mn In equation (B.7) Is a M x M matrix. The evaluation 
of R mn Is time consuming and It takes the most computer time In the analysis. 
Reference 10 shows a closed form solution to equation (B.7) for special case of 
the array element located In the x-y plane having Identical polarization 
parameters . 


DIRECTIVITY 


The directivity Is defined by 


D(e,(p) = 
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D(e,(p) Is known as total directivity. Also the reference directivity and 
the cross directivity can be easily obtained (ref. 14). 


Reference directivity: 


Cross directivity: 
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NUMERICAL RESULTS AND DISCUSSIONS 

This section presents some numerical results to demonstrate the applica- 
tions of the computer program. In order to substantiate the accuracy of the 
generalized array formulation and computer program, detailed comparisons were 
made with the results presented by King and Wong (ref. 2) and experimental 
data obtained at NASA Lewis (ref. 15). King and Wong reported on an N x N 
planar array configuration with symmetrical element patterns of the cos 
(ejd-type. They used direct Integration to compute the radiated power. 

The examples considered were an array of 2x2 elements for various element 
spaclngs and a 3x3 array for which element pattern and frequency were 
varied. In the NASA Lewis experimental case, a 2x2 array of rectangular 
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horns was used. In this case element pattern were measured In the array 
environment (to account for mutual coupling). 

Very good agreement was obtained In all cases. The 2x2 array 
reported by King and Wong assumed a symmetrical E-H-plane patterns with a q = 
3.54. Figure 4 shows a graphical description of the directivity as a function 
of element spacing for this array. The 3x3 array example used a sym- 
metrical pattern but the performance Is described as a function of element 
pattern, (cos (e)1, varying q) and operating frequency. Table II shows 
the results from both approaches. In the NASA Lewis experimental case the 
2x2 array was Investigated relative to far-fleld patterns. These patterns 
were measured at three different scan angles (boreslght, 3° and 5°). The ele- 
ment spacing and frequency were fixed (S = 2.5 frequency = 30 GHz). These 
results are presented In figures 4(a) to (e). 

A user guide for the programs developed Is presented In appendix A. The 
Implementation of equations (A1.3), (B.8), and (C.l) (antenna pattern, power, 
and directivity, respectively) with a computer program Is given In appendix B. 
This program (appendix B) can be easily Interfaced with available plotting 
routines for displaying the far-fleld antenna patterns. The numerical tech- 
nique used to solve for equation (B.8) Is not unique In any sense, but It was 
found to be faster than just using direct Integration. Many other techniques 
can also be used, and easily Implemented In the computer program (appendix B) . 


CONCLUDING REMARKS 

One of the advantages of this generalized array formulation Is that Is 
does not break-down for special cases as might occur In approximations using 
closed forms. Also the formulation developed does not limit the pattern obser- 
vation to the far-fleld zone. This can be very useful If the generalized for- 
mulation Is going to be used with analysis programs for dual reflector con- 
figurations. The program developed can be easily modified to be Implemented 
In the analysis of reflectors with phased array feeds. 

This computer program Is one of the key research tools at the NASA Lewis 
for analyzing advanced space communication antenna systems. The generalized 
formulation and computer program provides complete flexibility In analysis of 
array configurations and In the accurate analysis of experimental data. 
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APPENDIX A 


USER GUIDE 

PACALl 


Program Description 

Given an angle phi and an array of sources, each with its 
current magnitude and relative phase, this program 

calculates the cross polarization, reference polarization, 
and the far field magnitude for a series of angles theta. 

Input (FT05) 

X,Y, Z: coordinates of each source (in meters). 

AMPL: current amplitude for each source (in amperes). 
APHA: relative phase for each source (in radians). 

H: total number of sources. 

POL: denotes polarization 

1: linear polarization (x-direction) 

2: linear polarization (y-direction) 

3: right-hand circular polarization 
4: left-hand circular polarization 
QE,QH: exponent of the cosine function that is used to 
approximate the element pattern in the analytic form. 
INZT: interval between each elevation angle theta. 

I: l+phi, parameter to set for desired cut. 

Extraneous Input: (to be changed accordingly) 
line 2000 AUAVE=3E8/ (frequency) 

line 2500 ZETA=(- (total range of 

angles)+(j-l)*INZT)*PI/180 

Output : 

(FTOO) 

AZETA: value of theta at which the field values are 

taken (units degrees). 

ERAD: far field magnitude. 

(FT07) 

AZETA: (see FT08) 

AREF: reference polarization magnitude 
(FT06) 

AZETA: (see FT08) 

ACR: cross polarization magnitude 

Using the Program , 

Create an input file assigned to FT05 and output files 
assigned to FT06,FT07, and FT08. Running the program will 
fill FT06,FT07, and FT08 with data that can be used in the 
program ZPL0T1 to plot the appropriate graphs, while also 
printing the number of points calculated. 
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To run: 


DDEF FT05F001 , VS, INI 
DDEF FT06F001,VS,PAOUT3 
DDEF FT07F001,VS,PAOUT2 
DDEF FT08F001, VS,PA0UT1 
PACAL1 


Example of INI: 

CINPUT X=. 015, -.015, -.015, .015 Y=. 015, .015, -.015, -.015 

AMPL=1, 1,1,1 APHA=0, 0,0,0 M=4,QE=22,QH=16,INZT=. 1,1=1 SEND 


PAA1 


Program Description 

This program takes a series of sources with different 
amplitudes and phases, and determines the power via 
two methods. 

Input (FT09): 

tl: number of sources XX,YY,ZZ: x,y and z coordinates of 

source (meters) A: amplitude of source PHI: relative 

phase of source. 

IPOL: denotex polarity 

1: X-linear polarized 
2: Y-linear polarized 
3: circularly polarized 
QE,QH: exponents of cosine functions.. 

Output: 

(FT15) YPOWER: the power using direct integration method. 

This serves as an input to PADTV. 


PADTV 


Program Description 

This program requires the input filb to PAA1 as well as the 
output file and then calculates the peak directivity. 

Input: 

(FT09) Input to PAA1 
(FT15): Output to PAA1 

Output (FT06) 

DIR: peak directivity 
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ZPL0T1 


Program Description 

This program is designed to plot either the far-field, 
reference polarization, or cross polarization as calculated 
in PACAL1. 

Input : 

(FT06) 

X(I),Y(«J) (meters) 

These values are the output of the previous program 
PACAL1. If the plot of the far field is desired, the 
file which was assigned to FT08 in PACAL1 should now 
be assigned to unit 6. Similarly, if the plot of the 
reference or cross polarization is desired, the file 
assigned to FT07 or FTOG, respectively, should now be 
assigned to FT06. 

Extraneous Input: 

lines 1200,1250 

IVARS=NP=total number of points-1 
(This can be obtained from the output of PACAL1) 
lines 3500,3600 

VARS(4)=lower boundary of angles to be plotted. 
VARS(5)=upper boundary of angles to be plotted. 
Using the Program 

After running PACAL1, there will be data in the files 
assigned to FT08,FT07, and FT('6, which in this case will be 
denoted as A,B,and C. Assign either A,B,orC, to FT06 (after 
releasing (FT08,FT07, and FT06), depending on whether the 
far field, reference, or cross polarization plots are desired 
. Then, run the program with the appropriate plotting 
device (sidecar 4015) and the plot will appear. 

To run: 


RELEASE FT08 

RELEASE FT07 

RELEASE FT06 

DDEF FT06F001,VS,PAOU7 1 

ZPL0T1 

If this routine is executed after that shown in the PACAL1 
section, the far field will be plotted. 
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PROCDEF PAPLOT1 


This procdef runs the programs PACAL1 and 2PL0T1 in 
succession while defining the necessary input and output 
devices. 

As it is listed, file PA114 should contain the input to 
PACAL1. After assigning the devices FT08,FT07, and FT06 to 
files PA0UT1, PA0UT2, and PA0UT3, respectively, running 
PACAL1 will fill the respective file with the far field, 
reference polarization, and cross polarization magnitudes. 
After running PACAL1, setting the device FT06 to PA0UT1 will 
cause ZPL0T1 to plot the far field magnitude. 

The Procdef: 

ERASE PA0UT1 ERASE PA0UT2 ERASE PA0UT3 DDEF 
FT08F001, VS,PA0UT1 DDEF FT07F001,VS,PAOUT2 DDEF 
FT06F001, VS,PA0UT3 DDEF FT05F001,VS,PA114 PACAL1 
RELEASE FT08 RELEASE FT07 RELEASE FT06 RELEASE FT05 
GRAPH2D DDEF FT06F001 , VS,PA0UT1 ZPL0T1 




n 



APPENDIX B 


0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 

0002300 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0004900 

0005000 

0005100 

0005200 

0005300 

0005400 

UUUD3UU 

0005600 

0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006700 


PCAL1 


CXXXXXTHIS PROGARM WILL CALCULATE THE TRANVERSE ELECTRIC FIELDXXXXXXXXX 
CXXXXXAND THE REFERENCE/CROSS POLARIZATION COMPONENTS OF THE FIELDXXXXX 


CxxxxxdIMENSION STAMENTSxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL INZT 

DIMENSION XC100>,Y(100),AMPL'100),APHA<100) 


CxxxxxiNPUT DATA: XCI) ; X-COORDINAT E OF HORN I,Y(I> ; Y-COORDINATE OFXXXX 
CxxxxxHORN I.AMPL(I) ; AMPLITUDE OF I HORN, APHA(I) ; PHASE OF THE I xxx 
CxxxxxHORN.M ; NUMBER OF HORNS, POL ; 1 =X-POL , 2=Y-P0L , 3-RHCP, 4=LHCP . xxxx 
CXXXXXQH ; EXPONENT H PLANE; QE ; EXPONENT OG E-PLANExxxxxxxxxxxxxxxxxx 


NAMELIST/INPUT/X,Y,Z,AMPL,APHA,M,POL,QH,QE, INZT, I 

READC5, INPUT) 

INPH=1 

J=0 

AWAVE=3E8/10E9 
PI=4 . XATANC 1 . ) 

AK0=2.XPI/AWAVE 

CxxxxxsET THE ANGLES FOR THE PLOTxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

310 CONTINUE 
J=J + 1 

ZETA=(-88.+CJ-l)XINZT)xpi/180. 

APHI=(I-1)*INPHXPI/180. 


AZETA=ZETAX180 ./PI 
AAPHI=APHI*180./PI 
IFCAZETA.GT.88. )GO TC il 
IFCAAPHI .GT. 180 . )GO TO 999 

RAZP=0 . 

AIMZP=0 . 

DO 20 K=1,M 

CxxxxxsTART THE ARRAY FACTOR SUMATIONXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ANG=APHA(K)+X(K)XSIN(ZETA)XC0SCAPHI)+Y(K)XSIN(ZETA)XSIN(APHI) 
ANG=AKOXANG 

RAZP=RAZP+AMPL(K)XC0S(ANG) 

AIMZP=AIMZP+AMPL(K)XSIN(ANG) 

20 CONTINUE 

CxxxxxsET POLARIZATION PARAMETERSxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
I F( POL . EQ . 1 . ) GO TO 500 
I F( POL . EQ . 2 . ) GO TO 550 
IF(P0L.EQ.3.)G0 TO 600 
IFCPOL . EQ.4. ) GO TO 650 
GO TO 999 
500 Al=l. 

A2 = 0. 

51 = 0 . 

GO TO 800 
550 A1=0. 

A2 = l . 

SI = 0. 

GO TO 800 

600 A1=1./SQRT(2.) 

A2 = l . /SQRT ( 2 . ) 

SI=PI/2. 

GO TO 800 

650 Al=l ./SQRT (2 . ) 

A2 = l . /SQRT ( 2 . ) 

SI=-PI/2. 
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0006800 

0006900 

0007000 

0007100 

0007200 

0007300 

0007400 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008400 

0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010300 

0010400 

0010500 

0010600 

0010700 

0010800 

0010900 

0011000 

0011100 

0011200 

0011300 

0011400 

0011500 

0011600 

0011700 

0011800 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012500 

0012600 

0012700 

0012800 

0012900 

0013000 

0013100 


800 CONTINUE 


CxxxxxSTART THE ELEMENT PATTERN CALCULATIONxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
AUE=(COS(ZETA) )XXQE 
AUH=(COS(ZETA) )XXQH 
IF(AUE.LT.lE-3)AUE=lE-3 
IF(AUH.LT.lE-3)AUH=lE-3 
REEZL=A1XC0SCSI)XC0S(APHI)+A2XSIN(APHI) 

REEZL=AUEXREEZL 

AIMEZL=A1XSIN(SI)XC0S(APHI) 

AIMEZL=AUEXAIEMZL 

REEPH=A2XC0S(APHI)-A1XSIN(APHI)XC0S(5I) 

REEPH=AUHXREEPH 

AIMEPH=-1XA1XSIN(APHI)XSIN($I) 

AIMEPH=AUHXAIMEPH 


exxxxSTART THE TOTAL FIELD CALCULATION AT ZETA PHI OBSERVATIONXxxxxxxxx 

AREZE=RAZPXREEZL-AIMZPXAIMEZL 

AIMZE=AIMZPXREEZL+RAZPXAIMEZL 

AREPH=RAZPXREEPH-AIMZPXAIMEPH 

AIMPH=RAZPXAIMEPH+AIMZPXREEPH 

AEZ=AREZEXAREZE+AIMZEXAIMZE 

AEP = AREPHXAREPH + AII1PHXAIMPH 

ERAD=AEZ+AEP 

CXXXXXSTART THE CALCU1 ATION FOR CROSS POL AND THE REFERENCE POLXXXXXXXX 
AB 1 =A1 XCOS ( SI )XC03( APHI )+A2XSIN( APHI) 

AB2=-1. XA1XSINCSI )XC0S( APHI) 

AB3=-1 . xAlxcOSC SI) XSINC APHI) +A2XC0S( APHI) 

AB4=A1XSIN(SI )XSIN( APHI) 

AB5=A1XC0S ( SI )XSINC APHI )-A2XC0S( APHI) 

AB6=A1XSIN(SI)XSIN(APHT) 

AB7 =A1XC0S ( SI )XC0S( APHI )+A2XSIN( APHI) 

AB8=A1XSIN(SI)XC0S(APHI) ! 


CxxxSTART THE CROSS AND REFERENCE COMPUTATIONxxxxxxxxxxxxxxxxxxxxxxxxxx 
REER=(AREZExAB1-AIMZEXAB2)+(AREPHXAB3-AIMPHXAB4) 

AIMER=( AREZEXAB2+AIMZExAB 1) + ( AREPHX AB4 + AIMPHX AB3 ) 
RECR=(AREZEXAB5-AIMZEXAB6 )+(AREPHXAB7-AIMPHxAB8) 
AIMCR=(AREZEXAB6+AIMZEXAB5)+C AREPHX AB8+AIMPHXAB7) 
AREF=REERXREER+AIMERX AIMER 
ACR=RECRXRECR+AIMCRXAIMCR 
WRITE(8,400)AZETA,ERAD 
WRITE(7,400)AZETA,AREF 
WRITE(6,400)AZETA»ACR 
400 F0RMATC5X, E15. 5, 5X, E15. 5) 


GO TO 310 

11 WRIT E( 12 / 888 ) J 

888 F0RMAT(5X, 'TOTAL NUMBER OF POINTS',15) 


999 STOP 
END 
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0000100 DIMENSION XX( 1 0 0 ) , YY( I 0 0 ) 

0000200 DIMENSION AMAT ( 1 0 0 , 1 0 0 ) , YAMAT ( 1 0 0 , 1 0 0 ) 

0000300 DIMENSION GGXH ( 1 0 0 ) , GGXEC 1 0 * ) 

0000400 DOUBLE PRECISION B E( 1 0 0 , 1 0 0 ) , BH( 1 0 0 , 1 0 0 ) , CEC 10 0 , 1 0 0 ) , CH ( 1 0 0 , 1 00 ) 

0000500 DIMENSION CDE ( 1 0 0 , 1 0 0 ) , CDH( 1 0 0 , 1 0 0 ) , PHI C 1 0 0 ) , A ( 1 0 0 ) , THETA ( 1 0 0 , 1 0 0 ) 

0000600 DIMENSION AR ( 1 0 0 ) , A I ( 1 0 0 ) 

0000700 DOUBLE PRECISION QE , QH , P ( 1 0 0 , 1 0 0 ) , XK , XXE , XXH 

0000800 DOUBLE PRECISION B ESEC 1 0 0 , 100 ) , BESHC 1 0 0 , 1 00 ) 

0000900 DOUBLE PRECISION ARHO 

0001000 REAL LAMDA 

0001100 DOUBLE PRECISION YBEC 1 0 0 , 1 00 ) , YBHC 1 0 0 , 1 0 0 ) , YCEC 1 0 0 , 1 0 0 ) , YCHC 10 0 , 1 0 0 ) 

0001200 DOUBLE PRECISION ZA(100),ZB 

0001250 DOUBLE PRECISION XTHETA , DTHETA 

0001260 DOUBLE PRECISION AP1.AP2 

0001270 DOUBLE PRECISION AF11.AF12 

0001300 C 

0001500 CXXXXXTHIS ARE THE INPUTS TO POWER CALCULATIONXxxxxxxxxxxxxxxxxxxxxxxxx 

0001700 C NAMELIST/INPUT/M, XX.YY.ZZ, A, PHI, IPOL,QE,QH 

0001800 READ (9, INPUT) 

0001900 FREQ=30E9 

0002000 LAMDA=3. E8/FREQ 

0002100 PI=4 . XATANC 1 . ) 

0002200 XK=2XPI/LAMDA 

0002300 

0002400 

00 02500 

0002600 

0002700 Cxx**XTHIS SECTION WILL CALCULATE THE DISTANCE AND ANGLE M, Nxxxxxxxxxxx 
0002800 DO 40 JA=1,M 

0002900 DO 50 IA=1,M 

0003000 IFC(XX(IA).EQ.XX(JA)).AND.CYY(IA).EQ.YY(JA)))GO TO 55 

0003100 AL1=(XX(IA)-XXC JA) )X(XX( IA)-XX( JA) ) 

0 0 03200 AL2 = (YY(IA)-YY(JA) ) x ( YYC I A )-YY( JA ) ) 

0003300 AB1=SQRT ( AL1+AL2 ) 

0003400 AA1=(XX(IA)-XX(JA))/AB1 

0003500 AA2=(YY(IA)-YY(JA))/AB1 

0003600 

0003700 THETACIA, JA)=AA1XAA1-AA2XAA2 

0003800 PC I A , J A ) =AB1 

0003900 GO TO 89 

0004000 55 PCIA,JA)=0. 

0004100 THETACIA, JA)=1 . 

0004200 89 CONTINUE 

0004300 

0004400 

0004500 50 CONTINUE 
0004600 40 CONTINUE 

0004700 CxxxxxTHIS SECTION WILL COMPUTE THE BESSEL AND GAMMA FUNCTIONSXXXXXXXXX 

0004800 XXE=QE+.5 

0004900 XXH=QH+ . 5 

0005000 CALL GMMMACXXE, GXXE, IER) 

0005100 CALL GMMMA (XXH , GXXH , I ER ) 

0005110 AF1 1 =XXE+1 

0005120 AF12=XXH+1 

0005130 CALL GMMMA C AF1 1 , ARE1 , I ER ) 

0005140 CALL GMMMA ( AF12 , ARH1 , IER ) 1 

0005150 RATE=GXXE/ARE1 

0005160 RATH=GXXH/ARH1 

0005200 DO 60 JB=1,M 

0005300 DO 70 IB=1,M 

0005400 ARHO=XKXP(IB, JB) 

0005500 

0005600 

0005700 BESECIB, JB)=BESJPC ARHO, XXE) 

0005800 BESHC IB, JB)=BESJPC ARHO, XXH) 

0005900 70 CONTINUE 
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0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006700 

0006800 

0007000 

0007100 

0007150 

0007155 

0007180 

0007200 

0007300 

0007400 

0007500 

0007600 

0007610 

0007650 

0007655 

0007659 

0007660 

0007680 

0007681 

0007700 

0008000 

0008300 

0008600 

0008900 

0009315 

0009320 

0009400 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010250 

0010255 

0010300 

0010400 

0010500 

0010600 

0010700 

0010800 

0010900 

0011000 

0011100 

0011200 

0011300 

0011400 

0011500 

0011600 

0011700 

0011800 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012500 

0012600 

0012700 

0012800 

0012900 

0013000 


60 CONTINUE 
C 

CXXXXXTHIS SECTION WILL CALCULATE THE BE(I,J) AND BH(I,J) COEFICIENTSXX 


DO 80 JH=1,M 
DO 90 IH=1,M 
DTHETA=PI/200 
YBEC IH, JH)=0 
YBHCIH, JH ) =0 
YCECIH, JH)=0.0 
YCHCIH, JH)=0 . 0 
IFCIH.EQ. JH)GO TO 91 
DO 800 I L =1 , 1 0 0 
XTHETA=DTHETA*(IL-1) 

ZA(IL)=ABS(COS(XTHETA)) 

ZB=XK*P(IH, JH)XSIN(XTHETA) 
IF(COS(XTHETA).LT.1E-10)GO TO 800 
IF (ZB .LT. IE-10) GO TO 806 

FACT1=ZA(IL)xx(2*QE) 

FACT2=ZACIL)xx(2xQH) 

AP1 = 0 . 

AP2=2. 

ABE1=BESJP(ZB, API) 

ABE2=BESJP(ZB, AP2) 

YBEC IH , JH ) =YBE( IH, JH)+FACT1*ABE1*SIN(XTHETA)XDTHETA 

YBHCIH, JH)=YBH(IH, JHHFACT2XAB E1XS INC XTHETA)XDTHETA 
YCECIH, JH)=YCE(IH, JH)+FACT1XABE2XSIN(XTHETA)XDTHETA 
YCHCIH, JH)=YCH(IH, JH)+FACT2XABE2X5INCXTHETA)XDTHETA 
GO TO 800 

806 YBECIH, JH)=YBE(IH, JH)+FACT1X5INCXTHETA)XDTHETA 
YBHCIH, JH)=YBH(IH,JH)+FACT2XSIN( XT HETA1XDTH ETA 
800 CONTINUE 

IF (IH .EQ. JH) GO TO 91 
BBl=(2XX(QE-.5) )XGXXE/( (XKXPCIH, JH))KXXXE) 
BB2=(2XX(QH-.5) )XGXXH/C (XKXPCIH, JH) )XXXXH) 

BECIH, JH)=BB1XBESE(IH, JH) 

BHCIH, JH)=BB2XBESH(IH, JH) 

GO TO 90 

91 BE(IH,JH)=. 5/XXE 
BHCIH, JH)=.5/XXH 
YBECIH, JH)=.5XRATE 
YBHCIH, JH)=.5XRATH 


90 CONTINUE 
80 CONTINUE 
C 

CXXXXXTHIS SECTION WILL CALCULATE THE CHCI,J> AND CE(I,J) COEFICIENTSXX 
C 

DO 100 JC=1,M 
DO 110 IC=1,M 
CEC IC, JC) =0 
CH(IC,JC)=0 
DO 120 KK=1 , 12 
K=KK- 1 
FACT = 1 

DO 130 L=1,KK 
130 FACT = FACTXL 

XE=FLOAT(K)+QE+l .5 
XH=FL0AT(K)+QH+1.5 
CALL GLIMMA (XE , GXE , I ER ) 

CALL GMMMACXH, GXH, IER) 

IF (PCIC.JC) .EQ. 0) GO TO 225 
CCl=GXXE/( (4XXK)XFACTXGXE) 

CC2=(XKxp(lc, JC) )XX(2XK) 

CDECIC, JC)=CC1XCC2 
AK=FLOAT (K) 

IF ( AMODC AK , 2 . ) .NE. U) CDEC IC , JC ) =-CDE( IC , JC ) 

GO TO 226 


15 



0013100 

0013200 

0013300 

0013400 

0013500 

0013600 

0013700 

0013800 

0013900 

0014000 

0014100 

0014200 

0014300 

0014400 

0014500 

0014600 

0014700 

0014800 

0014900 

0015800 

0015900 

0015940 

0015950 

0016000 

0016100 

0016200 

0016300 

0016400 

0016500 

0016600 

0016700 

0016800 

0016850 

0016900 

0017000 

0017100 

0017200 

0017300 

0017400 

0017450 

0017500 

0017600 

0017700 

0017800 

0017900 

0018000 

0018050 

0018100 

0018200 

0018300 

0018400 

0018500 

0018600 

0018620 

0018650 

0018700 

0018800 

0018900 

0018940 

0018960 

0019000 

0019100 

0019800 

0019900 

0019950 

0019960 

0020000 

0020050 

0020100 


225 CDECIC,JC)=0 

226 CE(IC,JC)=CE(IC,JC)+.5*CDE(IC,JC) 

IF (PCIC,JC) .EQ. 0) GO TO 235 
DDl=GXXH/( ( 4**10 *FACT*GXH) 

DD2=(XK*P(IC, JC) )**(2*K) 

CDHCIC, JC)=DD1*DD2 

IF (AM0D(AK,2.) .NE. 0) CDH( IC , JC) =-CDH( IC , JC ) 
GO TO 236 

235 CDHCIC, JC)=0 

236 CHCIC, JC)=CH(IC, JC)+.5*CDH(IC, JC) 

120 CONTINUE 

CHCIC, JC)=CH(IC, JO-BHCIC, JC) 

CECIC, JC)=CE(IC, JO-BECIC, JC) 

110 CONTINUE 
100 CONTINUE 
C 

c 

c CALCULATING THE POWER 
C 

ZREAL=0 . 

ZIMAG=0. 

YZREAL=0 
YZIMAG=0 
DO 310 ID=1,M 
AR(ID)=A(ID)XCOS(PHI(ID) ) 

310 AI(ID)=A(ID)*SIN(PHI(ID)) 

IFCIPOL.EQ.DGO TO 150 
I F( IPOL . EQ . 2 ) GO TO 151 
GO TO 152 

150 DO 320 J=1,M 
DO 330 1=1, M 

AMAT<I,J)=BE(I,J)+BH(I,J)+THETA(I,J)x(CH(I,J)- 
YAMATCI, J)=YBE(I,J)+YBH(I,J)+THETAC V I,J)*(YCH(I 
330 CONTINUE 
320 CONTINUE 
GO TO 160 

151 DO 220 J=1,M 

DO 230 1 = 1 , n 

AMAT(I,J)=BECI,J)+BH(I,J)+THETACI,J)X(CECI,J)- 
YAMATCI, J)=YBE(I, J ) +YBH ( I , J ) +THETA C I , J ) * C YCEC I 
230 CONTINUE 
220 CONTINUE 
GO TO 160 

152 DO 420 J=1 , M 

DO 430 1=1, M 

AMAT (I,J)=BE(I,J)+BH(I,J) 

YAMATCI, J)=YBECI,J)+YBH(I,J) 

430 CONTINUE 
420 CONTINUE 
160 CONTINUE 

DO 520 J = 1 , M 
DO 530 1=1, M 

PREAL=(AR(I)*AR( J)+AI(I)*AIC J) )*AMAT(I, J) 
YPREAL=(ARCI)*ARC J)+AI<I)*AICJ))*YAMAT(I,J) 
YPIMAG=(ARCJ)*AI(I)-AR(I)*AI(J) )*YAMAT(I,J) 
PIMAG=(AR( J)*AI<I)-AR(I)*AI( J))*AMAT(I, J) 
ZREAL=ZREAL+PREAL 
ZIMAG=ZIMAG+PIMAG 
YZREAL=YZREAL+YPREAL ' 

YZIMAG=YZIMAG+YF’IMAG 
530 CONTINUE 
520 CONTINUE 

ZREAL=( 1/120 . ) *ZREAL 
ZIMAG=(1/120. )*ZIMAG 
YZREAL=(1/120. ) *YZREAL 
YZIMAG=( 1/120 . )*YZIMAG 
P0WER=SQRT(ZREAL**2+ZIMAGX*2) 

YP0WER = SQRT(YZREALxx2+YZiriAG**2) 

C 


CEC I , J ) ) 

, J)-YCECI,J)) 


CHC I , J ) ) 

, J)-YCHCI,J)) 
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0020200 
0020300 
0 0 2 0 4 0 0 
0020450 
0020500 
0020600 
0020700 
0020800 
0020900 
0021000 

00211 00 
0021200 


0021210 

0021220 

0021230 

0021300 

0021400 

0021500 

0021600 

0021700 

0021800 

0021900 

0022000 

0022100 

0022200 

0022300 

0022350 

0022355 

0022358 

0022359 

0022400 

0022500 

0022600 

0022700 

0022800 

0022900 

0023000 

0023100 

0023200 

0023300 

0023400 

0023500 

0023600 

0023700 

0023800 

0023900 

0024000 

0024100 

0024200 

0024300 

0024400 

0024500 

0024600 

0024650 

0024750 

0024800 

0024900 

0024950 

0025000 

0025100 

0025200 

0025300 


C POWER OF SINGLE ELEMENT 

SINGLE=(QE+QH+1)/(6Q*(2*QE+1)X(2XQH+1)) 

RATIO=POWER/SINGLE 

YRATIO=YPOWER/SINGLE 


WRITE(6,50Q)SINGLE 

500 F0RMAT(5X, ’POWER OF A SINGLE ELEMENT= ' , IX , E15 . 5 ) 

WRITE(6 ,501)RATIO,YRATIO 

501 F0RMATC///5X, 'CLOSED FORM RATIO= ' , IX, E15 . 5 , ’ DIRECT INT RATIO=', 

IX , 1E15 . 5 ) 

WRITE(6,502)POWER.YPOWER 

502 F0RMATC///5X, 'CLOSE FORM POWER= ' , IX ,E15 . 5 , ' DIRECT INT P0WER=',1X, 
E15.5) 

WRITE (15,600) YPOWER 
WRITE (16,600) POWER 
600 FORMAT (E15.5) 

C*****SPECIAL CASES FOR CHECKING RESULTS PREVIOUSLY CALCUL ATEDxxxxxxxxx 


C*xxx*LARGE SPACING CASE 6 LAMDA OR GREATERxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

POWER1=0. 

DO 550 1 = 1 , M 
P0WER1=P0WER1+A(I) 

550 CONTINUE 

P0WER1=SINGLE*P0WER1 
I F( IPOL . GE . 3 ) GO TO 810 
GO TO 811 

810 P0WER1=2.*P0WER1 

811 CONTINUE 
WRITE(6,503)POWER1 

503 F0RMAT(5X, 'LARGE SPACING POWER= ' , IX , E15 . 5 ) 


CxxxxxSYMMETRIC PATTERN QH=QE POWER CALCUL ATIONXxxxxxxxxxxxxxxxxxxxxxxx 
RECF=0 . 

REDI=0 . 

RIMCF=0 . 

RIMDI =0 . 

DO 560 J=1 ,M 
DO 561 1 = 1, M 

CA22 = (AR(I)xAR(J) + AI(I)*AI( '))*BE(I, J) 

CA33 = (AR(I)XAR(J)+AI(I)XAI^))XYBE(I, J) 

CA44=( AR(J)XAI(I)-AR(I)xAI( J) )XBE(I, J)- 
CA55=(AR(J)XAI(I)-AR(I)XAI(J) )XYBE(I, J) 

RECF=RECF+CA22 
REDI =REDI+CA33 
RIMCF=RIMCF+CA44 
RIMDI=RIMDI+CA55 
561 CONTINUE 
560 CONTINUE 

P0LPCF = SQRT(RECF*X2 + RIMCFXX2)X( 1/6 0 . ) 

P0LPDI=SQRT(REDIXX2 + RIMDIXX2)X( 1/60 . ) 

IF(IPOL.GE.3)GO TO 710 
GO TO 711 

710 P0LPCF=2.XP0LPCF 
P0LPDI=2. XPOLPDI 

711 CONTINUE 
WRITE(6,510)POLPCF,P0LPDI 

510 F0RMAT(///5X, 'POWER CLOSE FORM SYM= ' , IX, E15 . 5 , ' POWER D I SYM=',1X, 

E15.5) 

STOP 

END 
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PADTV 


0000100 

0000200 

0000400 

0000500 

0000550 

0000560 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001105 

0001110 

0001115 

0001120 

0001125 

0001130 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 

0002300 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0004900 

0005000 

0005100 

0005200 

0005300 

0005350 

0005360 

0005370 

0005380 

0005400 

0005450 

0005460 

0005470 

0005480 


DIMENSION THETA (10 0) , PHI ( 1 0 0 ), ETHETA ( 1 0 0 , 1 0 0 ), XX( 100 ) 
DIMENSION YY(IOO) ,AMP( 100) ,PHASE(100) 

DIMENSION AR(100,100),AI(100,100),AA(100,100) 
DIMENSION XRADI (100,100), EPHI ( 1 0 0 , 1 0 0 ) 

DIMENSION UE(100),UH(100),XIL(100,100) 

DIMENSION A ( 1 0 0 ) 

C 

C 

C READING VARIABLES 

C NAMELIST/INPUT/M, XX, YY.ZZ.QE.QH, A, PHI, IPOL 

READ (9, INPUT) 

DO 130 1=1, M 
AMP ( I ) =A ( I ) 

PHASE( I ) =PHI ( I ) 

130 CONTINUE 

READ (15,131) POWER 

131 FORMAT (E15.5) 

C 

PI =4 . #AT AN( 1 . ) 

FREQ=29 . 5E9 
XLAMDA=3.E8/FREQ 
XK=2*PI/XLAMDA 
Z0=377 
C 
C 

C TESTING POLARITY 
C 

IF ( IPOL .NE. 1) GO TO 102 
C=1 
B = 0 
PSI = 0 
GO TO 104 
C 

102 IF (IPOL .NE. 2) GO TO 103 
C = 0 

B = 1 
PSI = 0 
GO TO 104 
C 

103 C=1/SQRT(2.) 

B = 1/SQRT ( 2 . ) 

IF (IPOL .NE. 3) GO TO 105 
PSI=PI/2 
GO TO 104 
105 PSI=-PI/2 
C 
C 

C SOLVING EQUATIONS 
C 

104 DO 100 1=1,100 
DO 101 J=l,100 

C 

THETA(I)=(FL0AT(I-l)/99)KPI-PI/2 
PHI(J)=(FLOAT( J-l)/99)xpi-PI/2 
C 

UE(I)=(COS(THETA(I)))X*QE 
UH(I)=(COS(THETA(I) ) )**QH 
C 

AA1=UE(I)*UE(I) 

AA2=(C*C)*(C0S(PHI( J) )*COS(PHI( J)) ) 
AA3=(B*B)*(SIN(PHI( J ) ) XSIN ( PHI ( J ) ) ) 

AA4 = 2.*CXB*(C0S(PSI)XC0S(PHI( J) )*SIN(PHI( J)) ) 

ETHETA ( I, J)=AA1*(AA2+AA3+AA4) 

BB1 =UH (I)XUH(I) 

BB2=CXCX(SIN(PHI( J) )*SIN(PHI(J))) 
BB3=(B*B)*(C0S(PHI(J))*C0S(PHI(J))) 

BB4=AA4 

EPHI ( I , J )=BB1*(BB2+BB3-BB4 ) 
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0005500 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006610 

0006620 

0006800 

0006900 

0007000 

0007100 

0007200 

0007300 

0007350 

0007400 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008400 

0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009500 

0009600 

0009650 

0009651 

0009700 

0009710 

0009800 

0009900 

0010000 

0010100 


C 

XILCI, J)=ETHETA(I, J)+EPHI(I, J) 

ARC I , J ) =0 
AI ( I , J ) =0 
C 

DO 110 K = 1,M 

CC1=XKX(XX(K)XSIN(THETA(I))#C0SCPHI(J))) 
CC2=XKX(YY(K)XSIN(THETA(I) ) XS IN ( PHI C J ) ) ) 

EEXP = CC1 +CC2 

DREAL=AMP(K)xCOS(PHASE(K)+EEXP) 

AIMAG=AMP ( K ) XSIN (PHASEC K ) + EEXP ) 

ARC I , J ) =AR ( I , J )+DREAL 
AICI, J)=AICI, JJ+AIMAG 
110 CONTINUE 
C 

AACI, J)=AI(I, J)XAI(I, Jl+ARCI, J)XARCI, J) 
XRADICI, J)=AA(I, J)*XIL(I, J) 

101 CONTINUE 
100 CONTINUE 
C 
C 

C SORTING 
C 

DO 120 IA=1 » 100 
DO 121 J A = 2 , 1 0 0 

IF ( XRADI ( I A , 1 ) .GE. XRADI ( I A , J A ) ) GO TO 121 
XA=XRADI ( I A , 1 ) 

XRADI CIA, 1)=XRADICIA, JA+1) 

XRADI(IA,JA+1)=XA 

121 CONTINUE 
120 CONTINUE 

C 

DO 122 IB=2,100 

IF (XRADICI, 1) .GE. XRADI(IB,1)) GO TO 122 
XB=XRADI(IB,1) 

XRADI(IB,1)=XRADI(1,1) 

XRADICI, 1)=XB 

122 CONTINUE 
C 

WRITE(8,973)XRADI(1,1) 

973 FORMAT (5X, ' THIS IS THE PEAK VALUE’ , 3X, E15 . 5 ) 
DIR=4 . XPIXXRADI (1,1 l/CPOWERXZO) 

DIR =10 .XALOGIOCDIR) 

WRITE (6,123) DIR 

123 FORMAT ( 5X ,' DIRECTI VITY= ’, E15 . 5 ) 

STOP 

END 
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ZPL0T1 


0000100 

0000101 

0000103 

0000200 

0000300 

0000400 

0000401 

0000402 

0000500 

0000600 

0000700 

0000800 

0000900 

0000901 

0001000 

0001001 

0001100 

0001200 

0001250 

0001310 

0001410 

0001510 

0001610 

0001710 

0001810 

0001910 

0002000 

0002100 

0002200 

0002300 

0002400 

0002410 

0002420 

0002500 

0002600 

0002700 

0002705 

0002710 

0002720 

0002740 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0004900 

0005000 

0005001 

0005002 

0005004 

0005100 

0005200 


CxxxxxTHIS PROGRAM CAN BE USED TO PLOT THE ANTENNA 
CXXXXXFAR-FIELD PATTERN .( E-PLANE OR H-PLANE CUTS) 

CxxxxxBY R.J. ACOSTA 

DIMENSION XC5000) , Y( 50 0 0 ) , I VARS ( 2 0 ) , VARS ( 20 ) , RTNARRC 2 ) 

DIMENSION XTITLE(5)»YTITLE(5) 

LOGICALX1 IAXIS 
INTEGERX2 N1 

DATA XTITLE/'ELEV' , 'ATIO' , ’N AN ' , 'GLE ','DEG.V 
DATA YTITLE/'RELA' , ' TIVE' , ' AMP ' , • L ITU * , ' DE V 
CALL TITLEC4,20,15,XTITLE) 

CALL TITLE(3il7,15»YTITLE) 


IVARS ( 1 ) =2 
I VARS ( 2 ) =176 0 
NP=1760 

N1=NP 


DO 15 J=1,NP 
READ(6,300)X(J),Y(J) 

300 F0RMATC5X, E15.5,5X, E15.5) 

15 CONTINUE 

CALL SCLBAK( . FALSE. , N1 ,Y,RTNARR) 
VMAX=RTNARR ( 2 ) 


DO 16 J=1,NP 
Y(J)=Y( J)/VMAX 
IF(Y(J).LT.1E-8)Y(J)=1E~8 
Y(J)=10.XALOG10(Y(J)) 

16 CONTINUE 


VARS ( 1 ) =9 . 

VARS ( 2 ) =8 . 

VARS ( 3 ) =0 . 

VARS ( 4 ) =-90 . 

VARS ( 5 ) =90 . 

VARS(6)=6 . 

VARS ( 7 ) =2 . 

VARS (8 ) =-l . 

VARS ( 9 ) =0 . 

CALL XAXISC-l. ,~1. ,VARS) 
VARS (2 ) =9 . 

VARS C 3 ) =90 . 

VARS(5)=0. 

VARS ( 6 ) =8 . 

VARS ( 4 ) =-80 . 

CALL YAXISC-l. ,-l. ,VARS) 
CALL GPLOT(X,Y, IVARS) 
CALL DISPLA(l) 

READ( 9 , 9 93 )XYZ 
993 FORMAT ( 1A1 ) 

CALL TERM 
STOP 
END 
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TABLE I. - POLARIZATION PARAMETERS 


Polarl zatlon 
type 

a m 

b m 

'I'm 

Llnear-X 

1 

0 

0 

Llnear-Y 

0 

1 

0 

RHCP a 

1/ 2 

1/ 2 

0.5 ir 

LHCP b 

1/ 2 

1/ 2 

-0.5 ir 


a RHCP (right-hand circular 
polarized) . 

b LHCP (left-hand circular 
polarized) . 


TABLE II. - COMPARISON OF DIRECTIVITY RESULTS WITH THOSE OBTAINED BY 

KING AND WONG (ref. 2) 


Frequency, 

MHz 

S/X 

Element 

pattern 

HPBW, 

deg 

Pattern, 
cosQ (0) 
QE = <1H 

Klng-Wong 

measured 

directivity, 

dB 

\ 

Klng-Wong, 

dB 

calculated 

NASA 
Lewis, 
dB cal- 
culated 

450 

0.687 

139 

1.11 

17.1 

■m 

m 

500 

.763 

las 


17.9 


wttsm 

600 

.916 

89.4 

■m 

18.8 , 



700 

1 .068 

94.0 

0.91 

18.0 

17.4 

wstm 

800 

1 .220 

94.0 

0.91 

17.4 

16.9 

17.25 

















Figure 1. - Geometry of the generalized phased array. 



Figure 2. - A typical element coordinate system (X', Y'.Z'I and the reference 
coordinate system (X. Y.Z). 




Figure 3. - Geometrical picture ot the coordinate transformation 
between element coordinate system and the reference coor- 
dinate system. 



Figure 4 - Compared directivities for a 2x2 array as 
function of element spacing. 



Far-field magnitude, dB Far-field magnitude. dB 



(a) Boresight, reference polarization far-field antenna pattern. 
(E-plane cut). 

Figure 5. 


Experimental 



(b) Scan case <0^ - 3. 2°), reference polarization far-field antenna 
pattern. (E-plane cut). 


Figure 5. - Continued. 



Experimental 



Figure 5. - Concluded. 
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